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We develop a perturbative framework with which to discuss departures from ex- 
act Lorentz invariance and explore their potentially observable ramifications. Tiny non- 
invariant terms introduced into the standard model Lagrangian are assumed to be renor- 
malizable (dimension < 4), invariant under 577(3) <g> 577(2) ® U(l) gauge transformations, 
and rotationally and translationally invariant in a preferred frame. There are a total of 
46 independent CPT-even perturbations of this kind, all of which preserve anomaly can- 
cellation. They define the energy-momentum eigenstates and their maximal attainable 
velocities in the high-energy limit. The effects of these perturbations increase rapidly 
with energy in the preferred frame, more rapidly than those of CPT-odd perturbations. 
Our analysis of Lorentz-violating kinematics reveals several striking new phenomena that 
are relevant both to cosmic-ray physics (e.g., by undoing the GZK cutoff) and neutrino 
physics (e.g., by generating novel types of neutrino oscillations). These may lead to new 
and sensitive high-energy tests of special relativity. 
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1. Introduction 

Experimental tests of Lorentz invariance have become remarkably accurate. To give 
a quantitative measure of this accuracy, one imagines adding tiny Lorentz-violating terms 
to a conventional Lagrangian. Experiments can test Lorentz invariance by setting upper 
bounds to the coefficients of these terms. One common choice |jj is to alter the coefficient 
of the square of the magnetic field in the Lagrangian of quantum electrodynamics: 

B 2 -> (1 + e)B 2 . (1.1) 

Among other effects, this term causes the velocity of light c, given by c 2 = 1 + e, to differ 
from the maximum attainable velocity of a material body, which remains equal to one. 
(Shortly we shall consider more general Lorentz-violating perturbations.) 

The perturbation ( |1 . 1| ) breaks Lorentz invariance. It is translationally and rotationally 
invariant in the frame in which we are working ( "the preferred frame" ) but not in any other 
frame. If the preferred frame is the one in which the cosmic microwave background is 
isotropic, tiny and calculable anisotropies should appear in laboratory experiments. High- 
precision spectroscopic experiments that fail to find such anisotropies set the bound 
|1 - c 2 | = |e| < 6 x 10" 22 . 

In a note published last year ||] we pointed out that a better bound could be obtained 
from a very different sort of experiment if c < 1. In this case a charged particle traveling 
faster than light rapidly radiates photons until it is no longer superluminal. Thus no 
primary cosmic-ray proton can have energy greater than M/y/l — c 2 = M|e| -1 / 2 , where 
M is the proton mass. Because primary protons with energies up to 10 20 eV are seen, we 
set the bound 1 — c 2 < 10 -23 , almost an order of magnitude stronger than the atomic- 
physics bound. High accuracy is obtained from high energy rather than high precision. 
Moreover our bound requires no assumption about our velocity in the preferred frame. 

This effect, which we call vacuum Cerenkov radiation, is absent below a characteristic 
energy and turns on abruptly once that energy is reached. Such is not always the case, as 
the following example shows. 

Let \& denote a set of n complex scalar fields assembled into a column vector. If we 
assume invariance under the U(l) group \I/ — > e _lA \l/, the most general free Lagrangian is: 



£ = d^*Zd^ - **M 2 #, 
1 



(1.2) 



where Z and M 2 are positive Hermitian matrices. We can always linearly transform the 
fields to make Z the identity and M 2 diagonal, thus obtaining the standard theory of n 
decoupled free fields. Consider adding to the Lagrangian the Lorentz-violating term: 



where e is a Hermitian matrix. If e does not commute with M 2 , there is no way to disen- 
tangle the fields. The single-particle energy-momentum eigenstates go from eigenstates of 
M 2 at low momenta to eigenstates of e at high momenta. In contrast to vacuum Cerenkov 
radiation, this effect, which we call velocity mixing, turns on gradually. Gradual effects 
allow one to obtain high accuracy by combining moderately high energies with moderately 
high precision. 

A more striking gradual effect appears if this system is minimally coupled to electro- 
magnetism. In this case, a meson can decay to a less massive meson plus a photon at a 
rate growing with the cube of the energy. Analogous terms in the standard model can 
drive the otherwise-forbidden decay [i + —>■ e + + 7 and the 0-0 transition K + — > n + + 7. 

This example shows that what meant by "high energy" in this context depends very 
much on the details of the system under consideration. For simplicity, suppose n = 2 and 
let M 2 = diag (M 2 ,M|). The transition from eigenstates of M 2 to those of e occurs at 
energies ~ (M 2 — M 2 )/e\2- For the neutral kaon system, this energy is many orders 
of magnitude less than the characteristic energy of vacuum Cerenkov radiation (if the 
dimensionless invariance- violating parameters in the two processes are comparable). 

These are just illustrative examples. In §2 we study all local Lorentz-invariance vio- 
lating interactions that are rotationally and translationally invariant in a preferred frame 
and of renormalizable type (i.e., having mass dimension < 4).0 (Some of the results in §2.1 
and §2.2 were first established by Kostelecky and Colladay ||. We derive them anew here, 
both for completeness and because we wish to emphasize features relevant to high-energy 
tests.) 

Lorentz-violating perturbations can be divided into two classes, depending on whether 
they are even or odd under CPT. For a state with energy E, we show in §2.1 that the 

1 The condition of renormalizability can be given the usual justification: If we assume that 
breaking of Lorentz invariance occurs at some very small distance scale, only the renormaliz- 
able interactions survive the renormalization-group running of the couplings to experimentally 
accessible scales. 
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expectation values of the CPT-even interactions grow like E 2 for large E, while those of 
the CPT-odd interactions grow like E. Because we are interested in effects of very weak 
interactions made detectable by high energies, we limit ouselves primarily to a study of the 
the CPT-even interactions)! (Both our examples are of this class.) In §2.2 we construct 
the most general CPT-even interaction for the standard model and show that it preserves 
anomaly cancellation. In §2.3 we discuss certain approximations that we will use later on. 
As an illuminating exercise, in §2.4 we work out the kinematics of n-body decays in the 
special case in which off-diagonal matrix elements of the velocity-mixing matrices may be 
neglected. Novel phenomena arise. For example, a decay can be kinematically allowed 
both at low and high energy, but forbidden for an intermediate range of energies. 

The last section applies our formalism to various possibly observable manifestations of 
Lorentz violation. We discuss phenomena involving charged leptons in §3.1, in particular 
the possible appearance of radiative muon decay at high energies. We discuss phenomena 
involving neutrinos in §3.2, where we show how searches for neutrino oscillations at high en- 
ergy and long baseline can provide new and powerful tests of special relativity. We discuss 
hadronic manifestations of Lorentz-violation in §3.3, especially those relating to ultra-high 
energy cosmic rays. An Appendix explains why, contrary to our earlier assertion ||], the 
observed absence of a velocity difference between right- and left-handed photons (which 
would violate CPT) does not constrain the appearance of CPT- violating effects elsewhere. 

2. Generalities 

Here we develop the formalism needed to determine the observable consequences of 
Lorentz violation so as to obtain precise high-energy tests of special relativity. 

2.1. Building Lagrangians 

We first construct all CPT-even Lorentz- violating rotationally-invariant perturbations 
for a general renormalizable theory of scalars, spinors and gauge mesons. We also show 

2 This argument would be evaded if the CPT-even couplings were on the order of the squares 
of the CPT-odd ones, expressed as dimensionless ratios to an appropriate mass scale. This would 
be consistent with both renormalization-group flow and our energy-growth rule, and would lead 
to the dominance of the CPT-odd interactions at moderately high energies. We mainly ignore 
this possibility here. 
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that the matrix elements of these interactions grow with energy more rapidly than those 
of the CPT-odd Lorentz-violating interactions. 

We begin by summarizing some well-known properties of the Lorentz group SO (3, 1) 
|J. We assemble all the fields in the theory into a big vector $. The action of A, an 
element of 0( 3, 1), on these fields is effected by a unitary operator in Hilbert space U(A): 

U(A)^(x)U(A) = L>(A)$(A- 1 a;), (2.1) 

where D(A) is some finite-dimensional representation of the Lorentz group. 

The Lie algebra of SO(3, 1) may be written as the (complex) sum of two commuting 
angular momenta, J*-" 1 "-* and J^~\ An irreducible finite-dimensional representation of the 
group may be labeled by two half integers, j+ and j-, and is of dimension (2j + + l)(2j_ + 
1). We sometimes write the fields forming the basis for the irreducible representation as 
&m + , m_j where m± is the eigenvalue of Jz . 

The values of are (0,0) for a scalar; (0,1/2) or (1/2,0) for a Weyl spinor, 

depending on its chirality; (1/2, 1/2) for a 4-vector; (1, 1) for a traceless symmetric tensor; 
and the direct sum of (1, 0) and (0, 1) for an antisymmetric tensor. The complex conjugates 
of a set of fields transforming according to transform according to 

For R(e9), a rotation about an axis e by an angle 9, 

D{R(e9)) =exp [-i(J^ +) + J ( " } ) ■ e9]. (2.2) 

For B(e4>), a Lorentz boost in a direction e by rapidity (f>, 

D(B(ec/>)) = exp [(J (+ ) - J<H) . e>] . (2.3) 

Finally, for the anti-unitary CPT operator O, 

_1 $(x) O = (-l) 2j +&(-x). (2.4) 

We are ready to begin our analysis. From Eq. (|2.2|) we note that every rotationally 
invariant term £' in the Lagrangian must lie in a representation for which j+ = j- = j. 
Elementary angular-momentum theory tells us that the term is: 

j 

C! oc (-l) m $-,- m . (2.5) 

m=-j 
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If the state \ip) is boosted in the ^-direction by rapidity 0, the expectation value of £ is 
transformed according to: 

^\U\B{e z <ft) C'(0)U(B(e z <P))\^) oc e 2 ^(V# jW (0)|V>> + 0(e< 2 '"" 2 >*). (2.6) 

That is, (£') grows at large energy like -E 2 - 7 . As we shall show shortly, the largest value of 
j attainable with operators of renormalizable type is j = 1, a traceless symmetric tensor, 
which is CPT-even. The only other rotationally-invariant possibilities are j = 1/2, a 
vector, which is CPT-odd, and j = 0, a Lorentz invariant scalar. 

Let us begin by considering only scalar fields. With no loss of generality we can 
consider these to be all real. To attain j = 1 we need at least two derivative operators, and 
for renormalizability we can have no more than two (and no more than two scalar fields). 
Thus the only possibility is 

i J>(Te ab <9> & , (2.7) 

a, b 

where e a & is a real symmetric matrix and the sum runs over the scalar fields. (Of course, we 
could just as well have said that the only possibility is \ J2 a h do4> a e a bd (t) b ] the difference 
is a Lorentz invariant.) 

Now let us consider spinors. With no loss of generality we may take our fundamental 
fields to be a set of n (1/2, 0) Weyl spinors, u a ; the conjugate fields , « at , are then (0, 1/2) 
spinors. The only Lorentz-invariant interactions either couple two spinors of the same type 
or couple a spinor, a conjugate spinor, and a derivative. The most general free Lagrangian 
is 

uU(d - a ■ d)Zu + \u T o v Mu + \u ] a y M ] u*, (2.8) 

where Z is a Hermitian n x n matrix and M is a (possibly complex) symmetric matrix. 
We can always make a linear transformation such that Z is one and M is real, positive, 
and diagonal. We then have the theory of n free spin-1/2 particles, each of which is its 
own antiparticle. 

The only way to construct a renormalizable interaction with j = 1 is to couple together 
a spinor, a conjugate spinor, and a derivative. Thus we obtain 

±J2 te *bU at a-du\ (2.9) 

a, b 

where e a b is a Hermitian matrix. 
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A particularly simple case is that of two (1/2, 0) Weyl spinors carrying opposite charges 
under a (7(1) internal symmetry group. Then one Weyl spinor and the conjugate of the 
other can be joined to make a single Dirac bispinor, ip, the most general free Lagrangian 
is the standard Dirac Lagrangian, and the most general j ' = 1 interaction is 

\i,id ■ 7[e+(l + 75) + e_(l - 75 )]V, (2-10) 

where e± are real numbers. 

We now turn to gauge mesons. The couplings of the gauge mesons to the scalars and 
spinors is completely determined by gauge invariance, so we need only look at the gauge- 
field self-couplings. Let us begin with a single Abelian gauge field, and let us express it in 
terms of an electric and a magnetic field, defined, as usual, by E l — F° l and B % = |e y '*Fjft- 
Out of these we can make three independent rotational invariants of renormalizable type,S 
which we may choose to be E 2 — B 2 , E ■ B, and B 2 . The first two are Lorentz invariants; 
we are left with B 2 , just as in the first example of §1. 

A general gauge group is locally the product of simple factors and Abelian factors. 
Gauge invariance forbids cross terms between gauge fields belonging to different simple 
factors, so for each simple factor we have one interaction, of the form Yl a B a ■ B a , where 
the sum runs over the generators of the factor. For the Abelian factors we can have 
cross terms, but we can eliminate them by an orthogonal transformation on the Abelian 
generators; thus we again obtain one B 2 term for each factor. 

2.2. The Almost Lorentz-Invariant Standard Model 

As an example, let us construct the j • = 1 interactions in the standard model with 
three generations of quarks and leptons. There are only a few terms in the bosonic sector 
of the model: one B 2 term for each factor of the gauge group and one velocity-mixing term 
for the Higgs doublet, for a total of four real parameters. 

The number of parameters is much greater in the fermionic sector of the model. 
Each family of spinor fields transforms like the direct sum of five inequivalent irreducible 
representations of the SU (3) (g) SU (2) ®U(1) gauge group. Gauge invariance excludes cross 
terms between different irreducible representations, but not those between families. Thus 
we have five 3x3 Hermitian velocity-mixing matrices of the form fl2.9|) , for 45 additional 
real parameters. 

3 The CPT-odd term A ■ B, where A is the vector potential, is discussed in the Appendix. 
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A small reduction in the number of parameters may be effected by field redefinitions. 
We can rescale the space coordinates while leaving the time coordinate unchanged, thereby 
eliminating one of the gauge-field B 2 terms. Finally, we note that the minimal standard 
model conserves each lepton flavor. Thus, we may multiply the lepton fields by a phase 
factor depending only on flavor so as to make two of the off-diagonal elements in one of the 
lepton velocity-mixing matrices real while not affecting the Lorentz-invariant Lagrangian. 
Thus the standard model involves a total of 4 + 45 — 3 = 46 independent Lorentz- violating 
(but CPT-conserving) parameters. (We remark in passing that there are a comparable 
number of parameters corresponding to j = 1/2 CPT- violating departures from Lorentz 
invariance.) 

We used gauge invariance throughout this construction, but it has been the gauge 
invariance of the classical Lagrangian. Is this theory gauge-invariant as a quantum theory? 
Phrased another way, we know that anomalies cancel in the standard model, but does the 
cancellation persist when we take account of the 46 noninvariant terms? We demonstrate 
below that it does.B 

Our demonstration depends on the classic analysis of anomalies in Lorentz-invariant 
non-Abelian gauge theories 0. This begins by showing that the full theory is anomaly- 
free (to all orders in perturbation theory) if the corresponding theory of massless spinors 
coupled to external c-number gauge fields is anomaly-free. If we assemble all the spinor 
fields into a column vector u, the Lagrangian for such a theory is: 

ivt (D - a ■ D)u, (2.11) 

where the covariant derivative is defined by _D M = <9 M + A^, with a matrix- valued 
field composed of the gauge fields with their associated coupling matrices. If we define a 
gauge transformation of the fields in the usual way, 

5u = 5u)u, 5A lt = [5u,A ll ]-dp5w, (2.12) 

where 5uj is an infinitesimal gauge transformation, it may appear that the theory defined 



by ( 2.11 ) is gauge-invariant. However, this is not necessarily so. A carefully regulated 



computation of W(A), the generating functional of connected Green's functions, yields 



SW 



Tr / d 4 x e^8uF^F Xa , (2.13) 
3Z7v z J 



We were disappointed to discover this; we had hoped that the condition of anomaly cancel- 
lation would put some constraints on our 46 free parameters. 
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where F^ u = d^A v — d v A fl + L4 M , A v \. Only if this expression ("the anomaly") vanishes is 
the theory in fact gauge-invariant. Projection operators on irreducible representations of 
the gauge group commute with 5u and F^; thus a convenient way to evaluate the trace is 
to sum the contributions of the irreducible multiplets. For the standard model, this sum 
vanishes. 

We wish to extend all this to a theory with a Lorentz- violating interaction of the form 



|). That is to say, (|2.11| ) is replaced by 

iu*(D -(l-±e)ff-D)u, (2.14) 

where e is a Hermitian matrix acting only on the flavor indices. This replacement alters the 
high-momentum behavior of the spinor propagator and we must alter our regularization 
procedure accordingly. If we use massive regulator fields, the derivative terms in their 
Lagrangian must be of the form Q2.14] ) , not ( |2.11| ) . Likewise, if we regulate the measure in 



the path integral, we must use the differential operator from ( |2.14|) , not from (|2.11|) . 



We must sum over irreducible multiplets, which may be chosen to be eigenspaces of 
e. In each eigenspace, the Lagrangian is of the form ( |2.14j ) with e a number. We may 
introduce new variables (denoted by primes) by: 



X — ^ 1 q 6^ X ■, ^ — ' 1 1 



(2-15) 

A'={l-\e)A, A' Q = A . 

In terms of these variables, the e term disappears from the Lagrangian (and from the 
regularization procedure). The contribution of each multiplet to the anomaly is the same 
as it would be if e were zero, except that unprimed variables are replaced by primed 
ones. This replacement has no effect on ( |2.13p , which is invariant under general coordinate 
transformations. Thus the contributions of the irreducible multiplets are independent of 
e. If they cancel when e vanishes (as they do), they cancel for general e. 

2.3 From Lagrangians to Particles 

Particle properties of free fields follow trivially from the Lagrangian, but things are 
more difficult for interacting fields. For simplicity we study the energy-momentum relation 
for one real scalar field of (renormalized) mass m. The generalization to more complicated 
systems is straightforward. If the theory is Lorentz-invariant, the inverse renormalized 
propagator has the form: 

-iD' 1 = (p 2 - m 2 )A(p 2 ), (2.16) 



8 



for some function A. We normalize the field conventionally so that A{m 2 ) = 1, then add 
a Lorentz- violating interaction to the theory with some small coefficient e, as in Eq. ( |2.7|) . 
We begin in the linear approximation, retaining terms only first order in e. Later we 
investigate whether this approximation is justified. The addition to D~ x must transform 
like the 00 component of a traceless symmetric tensor. The only possibility is a multiple 
of 4p°p° — g 00 p 2 = Ap 2 + 3p 2 . The p 2 term can be absorbed in A, whence ( p.l 6|) becomes: 

-iD' 1 = (p 2 - m 2 )A(p 2 ) + ep 2 B(p 2 ), (2.17) 

for some function B. It is convenient to normalize the Lorentz-violating interaction such 
that B(m 2 ) = 1. 

To first order in e, the shift in the zeroes of D~ x is: 

p 2 = E 2 — p 2 = m 2 + ep 2 . (2.18) 

The energy-momentum relation may be rewritten in the seemingly conventional form E 2 = 
p 2 c 2 +m 2 c 4 , with c the maximal attainable velocity and mc 2 the rest energy of the particle. 
However, it must be remembered that c 2 has been changed by the factor to 1 + e and 
m 2 has been diminished by a factor of (1 + e) 2 . Of course, the tiny mass shift is of no 
experimental interest, but this is very much not the case for the shift in c 2 , as we saw in §1. 
Abrupt effects turn on when the dimensionless parameter ep 2 /m 2 is of order unity, while 
for gradual effects (like fi + — ► e + + 7) the energy at which the effect becomes significant 
can be many orders of magnitude smaller. 

Even for gradual effects, E is typically very large and we must ask whether new 
Lorentz-invariant physics at high energies might affect our predictions. Eq. (|2.18| ) shows 
that this does not happen. Even though E is large, p 2 remains 0(m 2 ) and the possible 
new physics remains irrelevant. 

When can the linear approximation be trusted? It can certainly be trusted for free 
field theory, where it is exact. It can also be trusted for leptons and electroweak gauge 
mesons. For these particles, all couplings are weak, all radiative corrections are small, and 
all propagators are well approximated by those of free field theory. 

Things are trickier for hadrons. A detailed investigation of QCD with Lorentz- 
noninvariant terms is beyond our ken, but we can make an educated guess on the basis of a 
simple model. Let's take QCD and rescale the space coordinates but not the time coordi- 
nate, as in Eq. ( |2.15|) . In the new coordinates, we seem to have a noninvariant theory, but 
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it's just QCD wearing a false beard. All that happens is that p 2 



in Eq. ( |2.16|) is replaced 



by p 2 — ep 2 . Thus an expansion in powers of e is in fact an expansion in powers of ep 2 /A 2 , 
where A is the QCD mass scale. Since typical hadron masses are 0(A), this implies that 
for gradual effects the linear approximation is a very good one, while for abrupt effects it 
is only a rough approximation. Of course, a rough approximation is not a useless one; it 
can give us a qualitative picture of what is going on, and even (with a modest amount of 
luck) yield correct order-of-magnitude quantitative predictions. 

2.4 The Kinematics of Particle Decays 

In this section we analyze the decay of a particle into n other particles in our Lorentz 
non-invariant theories ||. We make three simplifying assumptions: (1) that all particles 
are spinless. The extension to spin-| particles is straightforward and has no effect on our 
conclusions; (2) that the linear approximation is valid; and (3) that the matrix elements of 
the invariance- violating perturbation between particles with different masses are negligible. 
Thus we obtain a set of particles each of which obeys an energy-momentum relation of 
the form ([2.181) . That is to say, the ath particle has, in addition to its own mass, m a , its 
own maximum attainable velocity ("its own velocity of light") c a , and obeys the energy- 
momentum relation: 



In what follows we use a = for the decaying particle and a = 1 . . . n for the decay 
products. 

A decay is kinematically permitted if we can arrange the decay products such that 
their total momentum is po and their total energy is Eq. Let -E'min(Po) denote the minimum 
total energy of the decay products for given total momentum p . The decay is possible if 
and only if: 



because if E m i n < Eq we can obtain equality by adding opposite transverse components to 
two of the decay momenta. 

If we delete the transverse components of all decay momenta, we lower the final-state 
energy without changing the total momentum. It follows that all momenta are collinear in 
the configuration of minimum total energy. We use this fact to simplify our analysis and 
replace p a by p a , the longitudinal (and only nonzero) component of p a . 



tt2 ->2 2 , 2 4 




Eq > E m i n (po), 



(2.20) 
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E m [ n must be stationary under variations of the decay momenta that leave their sum 
unchanged. Introducing a Lagrange multiplier u, we must make 



J^Ea-ufy^Pa-Po) (2.21) 



stationary, where here and in the remainder of this section the sum is over the decay 
products. Differentiating with respect to p a , we find: 

u = ^=v a , (2.22) 

dp a 

where v a is the velocity of the ath particle. (We've used Hamilton's equations at the last 
step.) Thus all the decay particles move with a common velocity u. Furthermore, the 
relation: 

dE min = u^dpa = udpo , (2.23) 

shows that u = dE m i n /dpo. 

We can now explore the limits of small and large pq. For po = 0, the minimum energy 
configuration is one in which all the decay momenta vanish; the decay is allowed if and 
only if: 

m CQ > Y,m a c 2 a . (2.24) 

Of course, physics is certainly nearly Lorentz invariant, so that the c a 's can differ only 
very slightly. For all practical purposes, we can drop them from Eq. ( [2. 24]) . 

For very large po, u must be ultrarelativistic and we can approximate E a by c a p a . 
The energy is minimized by giving all the momentum to the particle with the smallest c. 
In this limit, the decay is allowed if: 

Co > minc a , (2.25) 

and is forbidden if Co < minc a . (We must go beyond the ultrarelativistic approximation 
to settle the question if cq = minc a .) 

Eqs. (2.24) and ( |2.25| ) are independent. A decay is allowed or forbidden in the high 



energy limit regardless of whether it is allowed or forbidden at low energies. It's interesting 
to ask what can happen at intermediate energies. For example, can there be alternating 
bands of allowed and forbidden energies? 
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We begin by answering this question when all the decay products have nonzero masses; 
afterwards we'll extend our analysis to deal with massless particles. Let's rewrite the 
condition for allowed decay, Eq > E m i n , as: 



m 2 4 > El in (p ) - c 2 p 2 = Y( Po ). (2.26) 



We shall prove that every stationary point of Y is a local maximum. Differentiating ( |2.2(j| ) 
yields: 

1 dY 

2dp~^ = EuiinU ~ c2 ° P ° = ~ ^ Pa ' < ^ 2 ' 27 ' 1 

where we have used E a = c 2 p a /u at the last step. The derivative vanishes at a stationary 
point, so at least one term in the sum must be negative, cq > min a ^o c a . Note that this 
condition coincides with Eq. ( |2.25| ), the condition that the process be allowed at high 
energy. 

Differentiating once more, we find: 

1 d 2 Y du v^/ 9 Osdp a du v^, 9 9n VaC 2 n 

At a stationary point, we can use the vanishing of ( [2.271) to write this as: 



1 d 2 Y du ^ 2 2 p a r 



C a C 



2 dp 2 dp ^ K a u/ U 

du ^ p a ( C 2 -C 2 ) 2 



|_ C 2 — U 2 Cq — W 2 



(2.29) 



u- 



E 



rfpo ^ (C 2 - W 2 )(cg - M 2 ) 



As po increases monotonically from zero to infinity, u increases monotonically from zero 
to minc a . Thus du/dpo is positive, as are the factors (c 2 — u 2 ) and (cq — m 2 ), these last by 
( |2.25|) . This completes the proof. 

Because every stationary point of Y is a local maximum, Y can have at most one 
stationary point. If Y has no stationary points, it is a monotone function of po and can 
cross TOq c o a t mos t once. This leads to four possibilities: (1) The decay is allowed at all 
energies. (2) The decay is forbidden at all energies. (3) The decay is allowed for all energies 
below a certain energy and forbidden for greater energies. (4) The decay is forbidden for all 
energies below a certain energy and allowed for all greater energies. If Y has one stationary 
point,i it may cross ttiqCq twice, the first time from below and the second time from above. 



5 The necessary and sufficient conditions for Y to have an extremum are Eq. fl2.25| ) and y^(c„- 
Cq) m a > (which ensures that Y(po) increases near po = 0). 
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This adds one more possibility: (5) The decay is forbidden for a certain band of energy 
but allowed for all energies above or below this band. There are no other possibilities. 

Our arguments break down if there are massless particles among the decay products. 
(For example, we can no longer write energy and momentum as functions of velocity.) 
Nevertheless, we now show that our conclusions remain valid. 

If there is more than one massless particle among the decay products, we can lower 
E m i n by giving all the momentum carried by the massless particles to the one with the 
smallest value of c. Thus no generality is lost by restricting ourselves to the case in which 
there is only one massless particle, which we label by a = 1. 

Consider a configuration with pi fixed and the remaining momentum, po — pi, dis- 
tributed among the massive particles so as to minimize their total energy. The total energy 
of this configuration is: 

E = c 1 p 1 +E , min (p -p 1 ), (2.30) 
where E' min is the minimum energy computed for the massive particles only. Thus, 

dE/dpi = c x -u'(p Q -pi), (2.31) 

where u' ', the common velocity of the massive particles, is a monotone increasing function 
of their total momentum. There are two cases: (1) u'(po) < c\, so that E is a monotone 
increasing function of p\ and E m i Q = E' min (po). The massless particle carries no momentum 
so that the analysis is the same as in the massive case. (2) u'(po) > ci, which is possible 
only if: 

ci < min c a , (2.32) 

for otherwise u' cannot reach c\. In this case, increasing pi lowers E until u' reaches ci, 
and 

Emi n = Ci(p "PO) + £min(Po)) ( 2 - 33 ) 

where p' is defined by u'(p' ) = c\. This analysis proceeds as in the massive case up to 
and including Eqs. (|2.27|) and ( |2.25|) , but the computation of d 2 Y jdp\ is different. Because 
only pi depends on po, 

\% = ^-® d t = { *-® Cl<0 ' <2 ' 34) 
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by Eqs. (|2.25|) and (|2.32|) . Just as in the case of only massive decay products, every 
stationary point of Y is a local maximum. Thus we reach the same conclusions here as we 
did there. i 



3. Applications 

Here we examine effects of CPT-conserving departures from Lorentz invariance on 
the behavior of ultra-relativistic particles. Our treatment of leptonic phenomena (in §3.1 
and §3.2) is not subject to the qualifications discussed in §1.4. However, the hadronic 
phenomena (in §3.3) often involve hadron energies so large that Ey/Sc 2 is comparable to 
the QCD mass scale. Thus, our treatment of the c a as energy-independent parameters 
is not always justified. This error affects the values of the Lorentz-violating parameters 
needed to produce novel phenomena, but not the phenomena themselves. 

3.1. Phenomena Involving Charged Leptons 

Tests of Lorentz Invariance from Photon Stability: Suppose that the velocity of light 
exceeds the maximal attainable velocity (MAV) of electrons, i.e., c 7 > c e , where for the 
moment we ignore possible flavor and helicity dependences. It follows [|J that photons of 
sufficient energy are unstable. In particular, the decay: 

7^e"+e + (3.1) 

becomes kinematically permitted when the photon energy E exceeds 2m e / v/5 7e , where 
5 7e = — Cg. The decay rate of photons well above this threshold is Y ee ~ ^ad ie E. The 
fact that primary cosmic-ray photons with energies up to 20 TeV have been detected lets 
us set the limit c 7 — c e < 10 -15 . Note that \S je \ is not so well constrained: the absence 
of vacuum Cerenkov radiation by electrons with energies up to 500 GeV sets the limit: 
c e — c 7 < 5 x 10 -13 . In a similar manner, much weaker constraints may be placed on the 
MAVs of the heavier charged leptons. 

6 The borderline situation, where the stationary point of Y occurs at po = p' , must be treated 
separately. In this case we must use the analysis for po < p'o to compute d 2 Y/dp^ from below 
and that for po > p' to compute it from above. The two answers don't agree because the second 
derivative of Y is not continuous at the stationary point, but it doesn't matter because both are 
negative. 
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More stringent tests of Lorentz invariance might be imagined to result from the sta- 
bility of photons under decay into two neutrinos, for which the threshold energy depends 
on tiny neutrino masses. This mode is forbidden in the standard model, which attributes 
neither masses nor magnetic moments to neutrinos, but there is a considerable body of em- 
pirical evidence for neutrino oscillations, and hence for neutrino masses. The non-standard 
physics responsible for neutrino masses could generate neutrino magnetic moments via loop 
diagrams, thereby enabling the decays 7 — > v + v' and v + v' (with v and v' necessarily 
distinct chiral neutrinos). 

The decay rate of a photon with energy E well above threshold into neutrinos is: 
^w' ~ ^ 2 ^b ! ) where the flavor-changing "magnetic moment" of the neutrino is \\xb 
with \x b = e/2m e (the Bohr magneton). We assume c v = c v > and put 5 1V = cl — c„. 
Because the magnetic interaction is dimension-five (rather than dimension- four, like the 
electric), T vv > is quadratic in 5 1V rather than linear .@ It follows that the range of an 
energetic photon is cosmological even if the neutrino magnetic moment is (implausibly) 
set equal to its experimental upper limit 0: A = 2 x 10~ 10 . Consequently, no strong 
bound on c 7 — c v can be deduced from observations of energetic cosmic-ray photons. 

Radiative Muon Decay? The decay mode mode \x — > e + 7 is often searched for 
but never found. This is not surprising: it is forbidden in the minimal standard model. 
Although induced by radiative corrections in models with neutrino masses, the expected 



branching ratio is far too small to be detected. However, and as we noted elsewhere |T(J [ 1 I 
Lorentz- violating perturbations that are not flavor diagonal can lift the accidental symme- 
try ordinarily preventing radiative decay. 

Recall that the velocity eigenstates of high-energy leptons do not in general coincide 
with their mass eigenstates at low energy. In the following analysis, we ignore possible 
mixings of electrons and muons with tau leptons. The relevant portion of the Lagrangian 
takes the following form in the preferred frame: 



U e)y.(p-el) {C fl i(l+ 75 )+C L i(l-75)} (ty 



(3.2) 



7 Our estimates of r ee and V vv i are lowest order in the appropriate Lorentz-invariant operators, 
with photons satisfying the Lorentz- violating dispersion relation E 2 — p 2 = 5E 2 ^"m 2 ". The 
decay rates are L ee | t ~ am 7 and V vv i | t ~ (A/is) 2 m 7 in the "photon rest frame." These 
expressions yield our results when boosted to the lab frame. 
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where \i and e denote fields corresponding to mass eigenstates. The matrices Cl,r (which 
would be unit matrices were Lorentz symmetry unbroken) are: 

r = I ( 2c a + 5c a cos 29 a Sc a sm29 a \ _ , , 

a "H 5c a sin2^ a 2c a -5c a cos2e a J ' a " ^ " ^ 

The mixing angle #l determines the velocity eigenstates of left-handed leptons (or right- 
handed antileptons) whose MAVs are cl ± ^Scl. Similarly, 9r determines the velocity 
eigenstates of right-handed leptons (or left-handed antileptons), whose MAVs differ by 
Scr. All four maximal velocities are known to be very close to c 7 , the velocity of light in 
vacuum. 

Electroweak gauge invariance implies that the matrix Cl appears in the kinetic-energy 
of neutrinos as well as charged leptons. In §3.2 we show how the parameters cl, Scl and 
9l may be constrained by experiments or observations involving neutrinos. 

For the analysis to follow, it is convenient to define the small parameters: 

e 2 = \5c a sin2# a | 2 , a = L, R. (3.4) 
The flavor-changing terms in Eq. (|3.2| ), namely 

| (cr + Q?) A* 7 • (P - eA ) e + \ -£r)U1 ■ (p~ eA) 75 e + h.c, (3.5) 

induce the decay process \i — * e + 7. We shall see that its helicity-dependent partial decay 
rate is ~ae^M7 3 , where M is the muon mass and 7 is its Lorentz factor in the preferred 
frame. This rate increases with the cube of the muon energy rather than falling with 
1/E. Thus flavor-changing Lorentz- violations, if they are present, will cause \x — * e + 7 to 
become the dominant decay mode of muons at sufficiently high energies. 

The rate of radiative muon decay is controlled by the muon energy and the magnitudes 
of the Lorentz-violating parameters. Otherwise, it is a first-order electromagnetic effect, 
not a weak decay. Its branching ratio is B ~ cte^7 Mtq, where tq is the lifetime of a muon 
at rest. Departures from Lorentz invariance also modify the rates of allowed processes, 
such as tv — > [i + v and \i — » e + v + P, but in these cases the conventional decay rates To/7 
and their Lorentz-violating corrections (of the form e 2 7 3 ro) are both intrinsically weak. 
They do not involve the enormous enhancement factor cxMtq ~ 2.6 x 10 15 that appears 
in the branching ratio for radiative muon decay. That's why the most sensitive tests of 
Lorentz invariance in this context are obtained from the study of muons, and in particular, 
from the search for a lifetime anomaly of muons at ultra-high energies. 
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The interaction ( j3.5| ), treated in lowest order perturbation theory, yields the rate for 
radiative muon decay. A straightforward but tedious computation (generously carried out 
for us by Mark Wise) yields the following result for its branching ratio when the muon is 
at rest: 

B = ^(el + el) * 6.4 x 10 14 (e| + 4) • (3-6) 



The current experimental limit [12 1 , B < 4.9 x 10 , yields an upper limit on the relevant 
Lorentz- violating parameter: 

e L + e R < 8 x 10~ 26 from muon decay at rest. (3-7) 

The branching ratio for Lorentz-violating radiative muon decay is a rapidly increas- 
ing function of the energy. Direct searches for it do not seem feasible. Nonetheless one 
might detect the onset of this mode through its effect on the muon lifetime, which for 
ultrarelativistic left-handed \i~ (or right-handed fx + ) is: 

rz,(7) = T ^ i , where b L = ^ (68 e\ + e 2 ) , (3.8) 

and for ultrarelativistic right-handed \i~ (or left-handed is: 

^) = TT^4> where ^ = ^(684 + 4)- (3-9) 

(These calculations were also carried out for us by Mark Wise.) At sufficiently high energy, 
the lifetime of muons with either helicity decreases with 7 -3 rather than increasing with 

7- 

The CERN g — 2 experiment, aside from measuring the muon's anomalous magnetic 
moment, offers a precise test of the energy dependence of its lifetime. At 7 = 29.3 (cor- 
responding to the 'magic energy' at which the experiment was performed) the results 
confirm the expected time-dilation of thed muon lifetime to an accuracy of one part in a 
thousand ||13|| . Because the muons in the ring favor neither helicity when time-averaged, 



we obtain the limit + < 2.7 x 10 , or: 

e| + t\ < 5 x 10" 25 from muon g-2. (3.10) 

which is inferior to that obtained from the direct search, Eq. (|3.7| ) , but not by much! (The 
agreement between theoretical and experimental values of g — 2 provides an independent 
but much weaker test of Lorentz invariance.) 
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The Muon Collider — A Threat Averted! Lorentz-violating effects consistent with the 
above constraints could interfere with the operation of a future muon collider. In [|10| and 



we took as a necessary condition for its proper functioning, that the muon decay 
rate at the design energy of the collider must not exceed twice its expected value. This 
criterion, assuming unpolarized beams, translates to \{bi, + &_r)7 4 < 1. Thus a 1 TeV 
collider requires b^ + &r < 2 x 10~ 16 (or e| + e 2 R < 3 x 10~ 32 ). The constraints discussed 
above are much weaker than this. However, we suggested that constraints sufficient to 
protect the muon collider might be obtained through studies of cosmic rays. 

The highest energy cosmic-ray muons arise from forward decays of secondary pions. 
These are mostly longitudinally polarised fi R and fi^, so that searches for a lifetime anom- 
aly of these muons can set a bound on 6r, which is proportional to the linear combination 
68 e i + e R- Although el can be severely constrained by observations of neutrino oscilla- 
tions, it's possible that e^. Absent any a priori knowledge of the ratio e^/e^, a 
bound on bn yields a 69-fold weaker bound on bu + bh 



In particular, we suggested |T0[ how a comparison of underground muon fluxes with 
different flight times but the same slant depth in rock might yield the required constraint. 
Following our suggestion, V.S. Narasimham and M.R. Krishnaswamy [[14] have analyzed 
data taken at the Kolar Gold fields. Their preliminary result is b^ < 8 x 10 -18 . An 
even stronger constraint has been obtained by R. Cowsik and B.V. Sreekantan |15| from 
their considerations of data concerning horizontal air showers. They obtain the remarkably 
strong bound 5r < 10 -25 . From Eq. (|3.9| ), this bound corresponds to: 

68 e 2 L +e 2 R < 10" 39 from air showers . (3.11) 

The results to either of these recent cosmic-ray analyses are sufficient to safeguard the 
muon collider from a weakening of the relativistic dilation of the muon lifetime. Indeed, 
taken by themselves, they provide sensitive tests of special relativity. For example, the 
air shower result (with = €r and maximal velocity mixing of muons and electrons) 
yields the following constraint on the difference of maximal velocities of the two velocity 
eigenstates: 

\c - c\ efi < 4 x 10" 21 . (3.12) 
We shall show in §3.2 how neutrino physics can do even better. 
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3.2 Phenomena Involving Neutrinos 

Although differences among neutrino velocities are severely constrained directly (by 
studies of neutrino oscillations, discussed below) and indirectly (from constraints involving 
muons or electrons discussed in §3.1 via 577(2) gauge invariance), there are only weak con- 
straints on neutrino maximal velocities as such. The current limit, \c v — c 7 | < 10 -8 , 



results from the detection of neutrinos from supernova 1987a |16|. A stronger limit, 
\c u — Cj\ < 2 x 10 -16 , may be set if neutrinos from gamma-ray bursts at cosmological 
distances could be detected Hl7j , but this result pales in comparison to other tests of 
Lorentz invariance that have been (or could be) set. 

Lorentz- Violating Neutrino Oscillations? We showed in || and [|IT| how CPT- 
conserving Lorentz violations lead to neutrino oscillations even if neutrinos are massless.i 
However, observable neutrino oscillations may result from a combination of effects involving 
neutrino masses and Lorentz violation. 

Neutrinos with modest energy (even solar neutrinos) are nevertheless ultra-relativistic 
particles because neutrino masses are known to be small. For this reason, searches for neu- 
trino oscillations can provide exquisitely sensitive tests of Lorentz invariance. We assume 
there exist three chiral neutrinos with Majorana masses given by the complex symmetric 
matrix m in a flavor-diagonal basis. Conventional neutrino oscillations are described in 
terms of the Hermitian squared-mass matrix m 2 = mrw . In particular, observable oscilla- 
tion effects depend on two differences of squared masses and four parameters akin to the 
Kobayashi-Maskawa angles and phase in the quark sector. 

With Lorentz symmetry violated, the description of neutrino oscillations becomes 
much more complicated. Although we usually neglect CPT-violating interactions, in this 
case it is hardly any work to include them. Thus, in addition to our usual velocity-mixing 
term, parameterized by c, a 3 x 3 Hermitian matrix of maximum attainable velocities, we 
allow the most general CPT-odd symmetry violating interaction, u^bu, where b is also a 
3x3 Hermitian matrix. The energies of ultra-relativistic neutrinos with definite momentum 
p are the eigenvalues of the matrix: 

cp + m 2 /2p + b. (3.13) 



Similar effects can result from violations of the equivalence principle rather than special 



relativity [18]. 
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Neutrino energy eigenstates in the limit of high energy are the eigenvectors of the matrix 
c, just as at low energy they are the eigenvectors of the matrix m 2 . 

To avoid undue complexity, we limit ourselves to a discussion of two-flavor neutrino 
oscillations. Imagine neutrinos to be produced with a definite momentum and flavor ( U£, 
where t = e or [i ) and detected after travelling a distance R through empty space. Their 
oscillations satisfy a seemingly conventional formula: 

P(u £ -»• vi) = 1 - sin 2 26 sin 2 {A.R/4} , (3.14) 

where the mixing angle and phase factor A appearing in Eq. (|3.14| ) (and expressed in 
the flavor basis) are given implicitly in terms of eight convention-independent parameters: 

A sin 26 = 5m 2 sin 26 m /E + 25b e iv sin 26 b + 25c e ir, 'Esm 26 c , 

(3.15) 

A cos 26 = 5m 2 cos 26 m /E + 25b cos 26 b + 25c E cos 26 c . 

The parameters characterizing the oscillations are: three mixing angles, 9 m , Q b and 9 C , two 
complex phases, r\ and rf, and the differences between the eigenvalues of the matrices m 2 , 
6, and c, denoted respectively by 5m 2 , 5b and 5c. (Note that 5c here is the same as 5cl in 
§3.1.) To illustrate the possibilities inherent in Eqs. ( |3.14j ) and ( |3.15| ), we mention a few 
special cases of Lorentz- violating two-flavor neutrino oscillations: 

sin \ (5m 2 R/4E) y/1 + (E/E ) 2 \ 

pi* -*)=i 1 1 + (W 1 <s-i«°) 



sin <^ (6v RE/2) y/l + (E /E) 2 

- 1 i + " (3 ' 1W) 

= 1 - sin 2 26 sin 2 {R(5m 2 /AE + 5b/2 + 5cE/2)} (3.16c) 



where Eq = 5m 2 /(25c). Eq. (|3.16a|) corresponds to 9 m = tx/A with 5b = sin26* c = 0. 
It yields maximal mass oscillations for E Eq, but essentially none for E ^> Eq. Eq. 
( |3.166|) corresponds to a converse case with 9 C = n/4 and 5b = sin2^ m = 0: maximal 
velocity oscillations at high energy, but none at low energy. 

To obtain Eq. (|3.16dj , we set all three mixing angles equal and put rj = rj' = 0. 



In this case, the energy-momentum eigenstates are independent of energy. This example 
encompasses all three scenarios discussed by Foot, Leung and Yasuda |L9| for atmospheric 
neutrino oscillations — each of which, they say, is consistent with the atmospheric neutrino 
data reported by Super-Kamiokande |2(| . 



20 



Conventional neutrino oscillations depend only on R/E, the ratio of the flight-length 
of the neutrino to its energy. If Lorentz symmetry is violated, the dependence on these 
parameters is more complicated. Nonetheless, neutrino experiments performed at a variety 
of energies can severely constrain the Lorentz-violating parameters. Let's give a simple 
example related to accelerator searches for v^-v e oscillations. The strongest current limit 
on Sm 2 (with sin 2 29 m ~ 1), 5m 2 < 0.09 eV 2 , follows from a relatively low energy experi- 
ment Higher energy neutrino experiments, such as [|22[| , offer less stringent constraints 
on 5m 2 but are better suited to search for Lorentz-violating velocity oscillations. From 
that experiment, and assuming sin 2 6 m ~ 1, we find for the difference of the maximal 
velocities of the two velocity eigenstates: 



c\ <6xl0 -22 . (3.17) 



Finally, we note that stringent constraints on the CPT-violating parameters that affect 
neutrino oscillations have been obtained from altogether different laboratory experiments. 



According to one of its collaborators ||23|| , the spectroscopic test of Lorentz invariance 
described in constrains the parameter 6 3 (as defined by Colladay and Kostalecky ||). 
They obtain | £>3 1 < 7 x 10 -19 eV for electrons and |&3| < 1.2 x 10 -21 eV for nucleons. The 
former result, expressed in our model and in the preferred frame, corresponds loosely! to 
the constraint |6 ee | < 3 x 10 -16 eV. This result suggests that CPT-violating effects are 
too small to affect neutrino oscillations, except when the source distance far exceeds the 
diameter of the Earth (as in the case of solar or extra-solar neutrinos). 

3.3 Phenomena Involving Hadrons 

To each particle species we assign a maximal attainable velocity c a . That is, we 
assume that a dispersion relation of the form E 2 = c 2 p 2 + m 2 c^ describes a particle of 
type a moving freely in the preferred frame. Many Lorentz-violating (but CPT-conserving) 
phenomena may be described in terms of the purely kinematic effects of these parameters. 
For simplicity, we ignore the helicity dependence of the MAVs, although it could easily be 
taken into account. Flavor-changing effects are not relevant to the phenomena discussed 
in this section and are likewise ignored. 

The Neutral Kaon System: Lorentz-violating effects can be abrupt or gradual. We 
gave examples of each in §1 — the sudden onset of vacuum Cerenkov radiation by energetic 



Here we neglect possible CPT violations involving right-handed electrons. 
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protons, and an energy-dependent modulation of the behavior of neutral kaons. Although 
our primary focus here is on abrupt hadronic effects, it is illustrative to examine the latter 
phenomenon in more detail. We consider the special case in which the velocity and mass 
eigenstates of neutral kaons coincide, and in which the MAVs of Kl and K$ are not 
the same: ck l — ck s 7^ 0. This leads to an energy dependence of their apparent mass 
difference, as determined by observations of time-dependent interference phenomena: 

AM = AM| rest + M^\c Kl -c Ka ), (3.18) 

where 7 and (3 are the Lorentz factor and velocity of the decaying kaons. A careful analysis 
of the experimental data carried out by Hambye, Mann and Sarkar |24| yields the bound:0 



\c Kl - c Ks \ < 3 x 10~ 2L . (3.19) 

Stable Neutral Pions? Let's turn to abrupt phenomena. A simple example involves 
neutral pions and photons. Suppose that c 7 — c^o > 0. The process it — > 27 (the 
dominant decay mode of neutral pions) becomes kinematically forbidden for pions with 



energies exceeding E = m^j -J — c*. Conversely, photons with energies significantly 
above E will decay rapidly according to the scheme: 7 — > 7 + 7r°. This example is not 
entirely academic. Suppose, for example: 

c 7 — c n o = c 7 — c e = c 7 — c^o = 10 -22 . (3.20) 

For this case, all modes of n° decay are kinematically forbidden at pion energies exceeding 
10 19 eV. Thus ultra-hig h-energy (UHE) primary cosmic rays may include neutral pions (if 
they are stabilized by tiny departures from Lorentz invariance), but not photons (which 
would be destabilized by the same mechanism). 

Stable Neutrons? Ordinarily, neutron decay (n — > p + e~ + v) is allowed but proton 
beta decay (p e + +n+v) is kinematically forbidden. As we have seen in § 2.4, departures 
from Lorentz invariance can affect the kinematics of decay processes. They even can invert 
this pattern! To see how this can come about, we examine the case c p = c e = c v < c n . 



10 Other tests of CPT and Lorentz invariance using neutral mesons are discussed in [24| and 
elsewhere |2{|. 
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Conventional relativistic kinematics may be used in this example (with c p as "the speed 
of light"), provided that the neutron is assigned an effective mass m e ff given by: 



2 _ 2 (2 2\ -2 



(3.21) 



where p is its momentum in the preferred frame. Neutron beta-decay is allowed if and only 
if m e ff > Trip + m e . Expressed in terms of the neutron energy E in the preferred frame, 
this condition becomes: 



E < Ei 



' mi — (Trip + m 



C p C n 



2.7 x 10 19 



10 



-24 



1/2 



eV. 



(3.22) 



With our choice of Lorentz-violating parameters, neutrons with energies exceeding E\ are 
stable particles that can be present among UHE cosmic rays. 

In a similar manner, we can deduce the necessary and sufficient condition for proton 
beta decay to be kinematically permitted. It is: 



E> E 2 ~ 



- (m p - m e y 



"mi 



4.1 x 10 19 



10 



-24 



1/2 



eV, 



(3.23) 



with E the proton energy in the preferred frame. For this partcular example, protons with 
energies exceeding E 2 are unstable particles that cannot be present among UHE cosmic 
rays. 

The above results are expressed in terms of a nominal choice, c p — c n = 10 -24 , lying 
beyond the sensitivity of current tests of Lorentz invariance. Thus it is conceivable that 
the highest energy cosmic-ray primaries are stable neutrons. 

Evading the GZK Cutoff? Soon after the discovery of the cosmic background radia- 
tion (CBR), Greisen []26| and Zatsepin and Kuz'min |27| saw how it limits the propagation 
of UHE cosmic rays. Primary nucleons with sufficient energy will suffer inelastic impacts 
with CBR photons. This results in what is known as the GZK cutoff, saying that nucleons 
with energies > 5 x 10 19 eV cannot reach us from further than ~ 50 Mpc. However, the 
primary cosmic-ray energy spectrum seems to extend well beyond 10 20 eV [28]. 

The mechanism producing UHE cosmic rays is unknown. Exotic origins have been 



suggested, such as: topological defects, active galactic nuclei, and gamma-ray bursts EJ. 
These schemes are constrained or ruled out by the GZK cutoff. Other explanations are 
designed to evade the GZK cutoff: a primary flux of magnetic monopoles [^0|, "Z-boson 
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bursts" produced by collisions of cosmic UHE neutrinos off relatively nearby relic neutri- 
lTTJ , and decay products of hypothetical super-heavy relic particles . 
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We have little to say about the origin of UHE cosmic rays. Rather, we point out that 
there may not be a GZK cutoff after all. Tiny departures from Lorentz invariance, too 
small to have been detected otherwise, have effects that increase rapidly with energy and 
can kinematically prevent cosmic-ray nucleons from undergoing inelastic collisions with 
CBR photons. The cutoff thereby undone, a deeply cosmological origin of UHE cosmic 



rays becomes tenable [33 



To see how the GZK cutoff is affected by Lorentz violation, consider the formation 
reaction yielding the first pion-nucleon resonance: 

p + 7(CBR) — >A(1232), (3.24) 

by which a nucleon with energy E collides inelastically with a CBR photon of energy uj. 
The target photon energies are a thermal distribution with temperature T = 2.73 K, or 
kT = ujq = 2.35 x 10 -4 eV. For a head-on impact, the condition Eq. (|2.20 ) determining 
whether the reaction can take place is approximately: 

Ml Ml 
2uJ+ 2E~ {CA ~ Cp)E+1 2E' (3 ' 25) 

where ca — c p is the relevant Lorentz-violating parameter. If ca = c p , Eq. ( |3.25| ) yields 
the usual threshold energy for this process, Ef = (M A — Mp)/4u. Otherwise it yields a 
quadratic inequality in E which can be satisfied if and only if (ca — c p ) < 5(u>) = u/2Ef. 
As ca — c p is increased toward <5, the threshold for the formation reaction grows toward 
2Ef. However, if it exceeds its critical value, 

2uj 2 

° A ~ ° p > Ml-M* ~ LT X 10 " 25 [U/UJ0]2 ' (3 ' 26) 



reaction (|3.24|) becomes kinematically forbidden for all E. Recalling that the photon 
spectrum is thermal, we see that if ca — c p ~ S(uq), the GZK cutoff due to resonant 
A(1232) formation would be relaxed. Should it much exceed this value, the formation 
reaction would be precluded off virtually all CBR photons. 

Reaction ( |3.24j ) is the dominant process leading to the GZK cutoff as originally for- 
mulated. However, if A(1232) formation is not possible, a weakened version of the cutoff 
may result from non-resonant photo-production: 

p + 7 (CBR) — > p + Ti . (3.27) 
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If ^ = c p ,, the threshold energy for pion production is E p = M 7T (2M p + M n )/4u. However, 
increased from zero, the threshold grows. As E — > oo, the pion energy .E^ 
must remain finite. Eq. (|2.20|) yields the kinematic condition for reaction (|3.27|) to occur: 



ml 



2oo > (c w -c p )E w + -^-. (3.28) 



This condition may be satisfied if and only if: 

2u 2 



c-x-CpK 6(u) e^~5x 10" 24 [u/ujo] 2 . (3.29) 



mi 



Reaction (|3.27 ), and multiple pion production as well, are kinematically forbidden at all 



proton energies if — c p > 5(u). For the actual case of a thermal gas of photons, 
c-K — Cp ~ S(uq) would suppress photo-pion production, or even eliminate it entirely so that 
no vestige of the cutoff survives. 

At present, lacking detailed observations of the highest energy cosmic rays and more 
precise tests of special relativity, we must regard as an intriguingly open question whether 
there is a GZK cutoff, and consequently, whether cosmic rays with energies above the 
cutoff can travel cosmological distances. 

Finally we note that tiny departures from Lorentz invariance, such as we have dis- 
cussed earlier in this Section, can explain the remarkable correlation discovered by Farrar 
and Biermann: that the five highest energy cosmic ray events appear to point toward 



compact radio- loud quasars |p4| . We suggest that these events could have been produced 
by UHE primary neutrons arising from sources at large redshift. These particles could be 
stable if c p > c n ; they could be immunized against the GZK cutoff if > c n . They are 
undeflected by intergalactic magnetic fields because they are neutral. 



4. Conclusions 

A wide variety of experiments and observations offer very precise tests of special rela- 
tivity. Many of these results can be interpreted in terms of differences of MAVs of different 
particles, such as would result from Lorentz- violating (but CPT-conserving) perturbations 

11 Note that much larger (and experimentally intolerable) violations of Lorentz invariance would 
be needed to produce a noticeable effect on the interactions of UHE cosmic rays with nuclei in 
the atmosphere. 



25 



of the standard-model Lagrangian. The strongest constraints of this kind of which we are 
aware have been mentioned earlier and are listed below: 



c p — c 7 < 1 x 10 

| Cm — C 7 | < 6 X 10 



-23 
-22 



<6 x 10" 22 



^3.2 ra 



\c-c\ <4xl0" 21 



\ck l -c k J <3 x 10 



-21 



j3.i m 

^3.3 m 



Two of these constraints result from cosmic-ray observations, the others from experiments 
performed at very high (accelerator) energies, or in one case, at low energy. They are 
consistent with strict Lorentz invariance to a precision of more than twenty-one decimal 
places. Possible symmetry violations, if present at all, must be exceedingly small. How 
much further must experimenters test special relativity; when is enough enough? Our 
analysis addresses this question. 

We have seen that maximal velocity differences lying two or three orders of magnitude 
below the current bounds can produce dramatic observable effects. Some would suppress 
or forbid the processes underlying the GZK cutoff, thereby permitting UHE cosmic-ray 
nucleons to travel cosmological distances. Another could stabilize UHE cosmic-ray neu- 
trons, which would point toward their astrophysical sources. Maximal velocity differences 
of neutrinos can help to explain the observed properties of both solar and atmospheric 
neutrinos. 

Fortunately, much can be done. Further observations of UHE cosmic rays are essential. 
They may confirm a predicted 'bump' just below the cutoff |35| resulting from products 
of inelastic collisions of primary protons with CBR photons, thereby providing evidence 
for the GZK cutoff. Or, they could belie the cutoff by confirming the Farrar-Biermann 



contention j3|J] that the highest energy events originate at cosmological distances. Dedi- 
cated searches for velocity oscillations of solar neutrinos, or of accelerator-produced ~ TeV 
neutrinos at baselines of ~ 1000 km, could reveal Lorentz-violating neutrino velocity differ- 
ences as small as 10 -25 . Finally, laboratory searches for diurnal anisotropies more sensitive 
than any done before have become feasible ||36||. 
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Appendix. Wrong Reasoning Made Right 

In our first paper on this subject Q we advanced an argument for neglecting CPT-odd 
Lorentz- violating interactions. 

The argument was based on studies [[57J of an electromagnetic interaction proportional 
to leoijfcA' 'F^ k = A - B. (This term is not gauge-invariant but it makes a gauge-invariant 
contribution to the action.) The new term makes an addition to the space-space part of 
the photon self-energy: 

Ely (p) -> Ely (p) + ive ijk p k , (A.l) 

where fi is a constant with dimensions of mass. The added term causes vacuum birefrin- 
gence; the experimental absence of this effect in radio-astronomy observations of distant 
quasars and radio galaxies leads to a very stringent upper bound, fx < 10~ 33 eV. In || 
we argued that any CPT-odd Lorentz-violating interaction would induce an addition to 
the Lagrangian proportional to A ■ B. Thus the extreme smallness of \x is evidence for 
the extreme smallness of all CPT-odd interactions, and further searches for CPT-violating 
effects are pointless. Unfortunately, this argument is invalid. 

Let us consider interactions for which the Lagrange density (not just the action) is 
invariant under electromagnetic gauge transformations; for brevity, we'll refer to these 
simply as "gauge-invariant interactions". We'll show that to first order in any gauge- 
invariant CPT-odd interaction, and to any order in the Lorentz-invariant interactions, 
\x vanishes. (This confirms the expectation that gauge-invariant terms in the Lagrange 
density cannot induce gauge- noninvariant ones.) Thus the fact that \x is known to be tiny 
offers no strong constraint on the magnitude of gauge-invariant CPT-violating interactions. 

It's useful to analyze a simple example before giving our proof in full generality. 
Let E Ml/ (p, q) be the one-particle irreducible (1PI) Green's function for two photons and 
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one Lorentz-violating interaction Lagrangian. Here the photon with index \i (y) carries 
momentum p (q). (Thus the Lorentz-violating interaction carries momentum — (p + q).) 
In first order in the Lorentz-violating interaction 



rv( P ) = iv(p,-p). (A.2) 

Because the Lorentz-violating interaction is assumed to be gauge invariant, F obeys the 
Ward identities 

P^vAPiQ) = 0, (A3a) 

and 

q v T liV (p 7 q) = 0. {AM) 

Now let us assume we compute V by summing only those diagrams that have no 
internal photons. These diagrams have only massive-particldll internal lines; thus, in this 
approximation, V is analytic at p = q = 0. Differentiating (A. 3a) with respect to p^ and 
setting p = 0, we find T^ u (0, q) = 0. Thus every nonvanishing term in the Taylor expansion 
of T has at least one factor of p; V is 0(p). The same reasoning applied to (A. 3b) tells us 
that T is 0(q). Since p and q are independent variables, V is 0(pq). Thus, by (A.2), the 
addition to the self-energy is 0(p 2 ) and makes a vanishing contribution to [i. 

The proof we have given rested on analyticity at vanishing momentum. Unfortunately, 
internal photons can produce singularities at precisely this point. In the remainder of this 
appendix we study these singularities and prove the following: 

Theorem: Let T^-.^niPi ■ ■ -Pn) denote the n-photon 1PI Green's function, and let us 
introduce a uniform scaling parameter A by replacing p a by Xp a . Then to either zeroth 
or first order in the CPT-odd interaction, and to any finite order in the Lorentz-invariant 
interactions, and for any positive e, vanishes as A goes to zero more rapidly than A n_e . 
(Note that our announced result, \x = 0, is a corollary of this theorem for n = 2.) 

The proof proceeds by induction in the number of internal photon lines. 

We begin with the case of no internal photon lines. In this case r^ n ) is analytic at 
vanishing external momentum, and, to first order in the Lorentz-violating interactions, 
the argument is a straightforward generalization of that for the simple example above. 
We treat the n external momenta as independent variables, letting the Lorentz-violating 

12 Neutrino masses may be tiny on the scale of high-energy physics, but they're enormous on 
the scale of radio waves. 
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interaction carry off the momentum inserted at the photon lines. We then use the Ward 
identities to show that the leading term in the Taylor series is 0(p\ . . .p n ) = 0(X n ). 

This doesn't work in zeroth order in Lorentz violation, because there is no Lorentz- 
violating interaction to carry off the momentum. The best we can do is to choose q± . . . q n -i 
as our independent variables and use the Ward identities to show that the leading term 
is 0(pi . . .p n -i) = 0(X n ~ 1 ). However, the term zeroth order in Lorentz violation is 
Lorentz invariant, and there is no Lorentz-invariant way to construct a rank n tensor as 
a multilinear function of n — 1 independent vectors. Thus this term must vanish, and the 
leading term is at best 0(X n ). 

To prove the inductive step, we need some information from the once well-known 



theory of Feynman-diagram singularities ||38|| . For real external momenta, the case of 
interest here, the theory can be reduced to a set of simple algorithms: (1) A given Feynman 
graph generates a family of reduced graphs, each obtained by taking some proper subset 
of the internal lines of the graph and contracting them to points. (2) The reduced graph 
leads to a singularity if it can be interpreted as a diagram of a classical process occur ing 
in space-time, with all particles (that is to say, all uncontracted internal lines) on the mass 
shell and moving forward in time. (3) If there is a singularity, the associated discontinuity 
is calculated by the usual Feynman rules except that (p 2 — m 2 + ze) _1 is replaced by 
5(p 2 — m 2 )9(p°) in the propagators for the uncontracted internal lines. 

We are now ready to use induction. We assume the theorem is true for r or fewer 
internal photon lines and consider the singularities of graphs with r + 1 internal photons. 
Since all energies inserted into the graph are arbitrarily small, there is not enough energy 
to create a massive particle and all massive-particle lines must be contracted. Thus the 
reduced graph contains only photon lines, joined together at vertices which represent con- 
tracted subgraphs of the original graph. If we sum, in some fixed order of perturbation 
theory, all graphs leading to the same reduced graph, the vertices become the T's at the 
appropriate order of pertubation theory. (We're being a bit sloppy here: subgraphs of a 
1PI graph need not be 1PI, so the vertices can contain tree graphs. However, it's easy to 
check that these have no effect on the power-counting of the next paragraph.) 

Because some of the internal photons in the original graph must be uncontracted, an 
m-line vertex in the reduced graph must have r or fewer internal photon lines and, by the 
inductive hypothesis, must vanish at small A more rapidly than A m_e . It will be convenient 
to consider this as one factor of A 1_e for each photon line attached to the vertex. We can 
now compute the A-dependence of the discontinuity. Every external photon contributes 
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a factor of A 1_e . Every internal photon contributes a factor of A 2_e from its two ends 
and a factor of A -2 from the 5-function, yielding no net contribution. Every independent 
loop integration contributes a factor of A 4 . Thus the discontinuity vanishes faster than 
A" +4L_e , where L is the number of loops. This in turn vanishes faster than A n_e , which 
is the result we need. (It is critical that we are computing the discontinuity and not the 
full T. The 5(p 2 )9(p°) propagators in the reduced graph keep the internal momenta small 
when the external momenta are small and legitimize the use of a small-momentum bound 
for the vertices.) 

Once we have the discontinuity, we can construct a function with that discontinuity, 
for example, by integrating a dispersion relation. This function also vanishes faster than 
A n " e . Thus r^ n ) is the sum of a singular function that vanishes faster than A n e and a 
function that is free of singularities, that is to say, an analytic function. But we can use the 
same arguments for the analytic function that we used for the case of no internal photon 
lines to show that it vanishes like A n . This completes the proof. 
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